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Abstract —In this paper, frequency-domain analysis based on 
frequency sweep method is presented for a nonlinear double 
integrator and a new linear integrator. All the two types of 
integrators can estimate the onefold and double integrals of a 
signal synchronously. With respect to the linear double integrator, 
the nonlinear integrator has better estimation performance and 
stronger robustness. Importantly, the Integrator parameters can 
be regulated from the frequency-domain analysis. 

Index Terms —Double Integrator, frequency-domain analysis, 
frequency sweep. 

I. Introduction 

This paper focuses on the frequency-domain analysis of two 
types of integrators, which can estimate the onefold and double 
integrals of a signal synchronously. 

Integrals are important components in almost all engineer¬ 
ing applications. The problem of double integral is that of 
estimating the number I 2 (a) = Jq/qO- (a) dadr with [ 0 , f] a 
finite time interval. Obtaining the double integral of a signal 
is crucial for many kinds of systems, especially for Inertial 
Navigation System (INS). 

The usual observers or differentiators [l]-[4] can estimate 
the derivatives of the signal, but not its multiple integrals. 
There are several linear approximated methods to estimate 
onefold integral [5-8]; Romberg integration, Gaussian quadra¬ 
ture, extended Simpson’s rule, fractional-order integrator. In all 
of the aforementioned studies, there is no stability analysis. 
Furthermore, they are easily disturbed by stochastic noise 
(especially non-white noise), and the drift phenomena occur in 
such systems. In [9], a fractional-order integrator is proposed 
to approximate the irrational fractional-order integrator I/s'". 
However, the condition of 0 < m. < 1 limits the application 
of the fractional-order integrator. Recently years, Kalman filter 
is used to handle the separation of probabilistic noise and to 
estimate the position and velocity from the acceleration mea¬ 
surement [10]-[12]. However, for Kalman filter, it is assumed 
that the process noise covariance and measurement noise 
covariance are zero-mean Gaussian distributed, and the process 
noise covariance is uncorrelated to the estimation error. These 
assumptions are different from the real noise in signal. The 
inaccurate noise information in sensed accelerations may lead 
to the estimate drifts of position and velocity. 

In [13], a nonlinear double integrator was presented based 
on finite-time stability [14, 15]. The proposed double inte¬ 
grator can estimate the onefold and double integrals of a 
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signal synchronously, and the stability and robustness in time 
domain were analyzed. The merits of the presented double 
integrator include its finite-time stability, ease of parameters 
selection, sufficient stochastic noise rejection and almost no 
drift phenomenon. The theoretical results are confirmed by 
an experiment on a quadrotor aircraft to estimate the position 
and velocity from the acceleration measurement. The nonlinear 
double integrator leads to perform rejection of low-level persis¬ 
tent disturbances. However, no robustness for high-frequency 
noise is analyzed, i.e., no frequency-domain analysis on the 
effect of high-frequency noise is considered. In fact, stochastic 
high-frequency noise exist in almost all signals. Therefore, 
with respect to the ability of rejecting of low-level persistent 
disturbances, the strong robustness of reducing high-frequency 
noise is also necessary. And the parameters selection based on 
frequency-domain analysis is required. 

In this paper, based on singular perturbation technique 
[16, 17], a linear double integrator is presented to estimate 
the onefold and double integrals of a signal synchronously. 
Moreover, using frequency sweep method, frequency-domain 
analysis is presented for the nonlinear double integrator in [13] 
and the linear integrator. From the frequency-domain analy¬ 
sis, comparing to the linear double integrator, the nonlinear 
integrator has the better estimation performance and stronger 
robustness. Also the observer parameters are more easily to 
be chosen from the frequency-domain analysis with respect to 
the analysis in time domain. 

H. Nonlinear and linear double integrators 
A. Nonlinear double integrator 

A nonlinear double integrator has been presented in [13], 
i.e., for system 

±1 = X2-,X2=X3; 

2 

e^xz = — kj [e'xil"* sign (xj) 

i=l 

-kslxs - a{t)\“'^ sign{x 3 - a(t)) ( 1 ) 

if input signal a (t) is the continuous and first-order deriv¬ 
able, then there exist 7 > 1 and F > 0, such that, for 
f > eF (S(e)e (0)), 

Icc,-a,(f)|<L£“'T'-',i=l,2,3 (2) 

where ai (t) = a (cti) daida 2 , 02 (t) = a (cti) dai; 

xi (0) = ai (0), X 2 (0) = 02 (0), X 3 (0) = 03 (0); £ e (0,1) 

is the perturbation parameter; L is some positive constant; 
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7 = (1 - /3)//3 > 1, and /3 e (0,ai/(ai +4)); ai,a 2 ,a 3 
satisfy: 


Q!3 € (0,1), q;2 


0:3 _ as 

2 — a 3 3 — 2a3 


ki,k 2 ,k 3 > 0 are selected such that 


(3) 


fci > 0,A:3 > 0,fc2 (4) 

kz 

ei = Xi -ai{t), i = 1,2,3; e =[ ei 62 63 E;(e) = 

diag{e, £^, £^}. 

In nonlinear double integrator (1), 0:3 tracks the input signal 
a{t), X 2 and xi estimate the onefold and double integrals of 
signal a{t), respectively. From Theorem 1 in [13], nonlinear 
double integrator ( 1 ) leads to perform rejection of low-level 
persistent disturbances. However, no frequency-domain analy¬ 
sis on the effect of high-frequency noise is considered. In fact, 
high-frequency noise exist in almost all signals. Therefore, the 
frequency-domain analysis for the nonlinear double integrator 
is inevitable. 

As we know, a linear system is easy to perform frequency- 
domain analysis with respect to nonlinear one. In the fol¬ 
lowing, based on the design of nonlinear double integrator, 
a simple linear double integrator will be designed (when 
as = 1), and Theorem 1 is presented as follow. 


B. Linear double integrator 
Theorem 1: For system 


Xi = X2-,X2=Xz] 

e'^xz = —kiexi — k 2 S^X 2 — kz {xz — a {t)) (5) 

if input signal a (t) is continuous, integrable and first-order 
derivable, then 

\m\Xi = tti (t) ,i = 1, 2,3 (6) 

E-S-O 

where oi (t) = f* fj a (s) dsdr, 02 (t) = f* a (t), 03 (f) = 
a (f); Xi ( 0 ) = Ui ( 0 ), i = 1 , 2 ; £ G ( 0 , 1 ) is the perturbation 
parameter; ki,k 2 ,kz >0 are selected such that 

ki > 0,kz> 0, k 2 > £^^ (7) 

kz 

In linear double integrator (5), xz tracks the input signal 
a(t), X 2 and a;i estimate the onefold and double integrals of 
signal a(t), respectively. 

Proof of Theorem 1: The Laplace transformation of Eq. (5) 
can be obtained as follow: 


sXi{s) = X 2 {s)-sX 2 {s) = Xz{s)-, 

2 

se^Xzis) = -J2kie^X,{s) 

-kz{Xz{s)-A{s)) ( 8 ) 

where (s) and A (s) denote the Laplace transformations 
of Xi and a (t), respectively, and s denotes Laplace operator. 
From ( 8 ), we obtain 


, X X, (s) 

^*(s) = ^^A = l,2,3,j€{l,2,3} 
Therefore, Eq. ( 8 ) can be written as 


(9) 




^ 2 ( 5 ) = 






sJ-^ 


Xj is) 


Then, it follows that 

^2 (s) 


A ( 5 ) 


-fc3M^-A(s) ( 10 ) 


( 11 ) 


s3-^+i£4 + + m 

2=1 


i.e., 


X, (s) 


^kz 


A is) 

Therefore, we obtain 


( 12 ) 


S^£4 -f X) 4fcj£* -I- s'^kz 
i=l 


y Xj (s) ,_3 

hm f, / = s-’ 4 


(13) 


£—^0 A is) 

where j G {1,2,3}. It means that Xi approximates Ui it) for 
1 < 2 < 3. 

Eurthermore, the denominator of Equation (12) is required 
to be Hurwitz, i.e., polynomial s _l_ is 

Hurwitz. It is equivalent that -f + k 2 S + ki should 
be Hurwitz. Eor arbitrary £ G (0,1), from the Routh-Hurwitz 
Stability Criterion, polynomial s^ + !^s‘^ + k 2 S + ki is Hurwitz 
if fci > 0,kz > 0, ^2 > s^ki/kz- This concludes the proof. ■ 

Erom Theorem 1, the presented linear double integrator 
can’t guarantee to perform rejection of high-frequency noise. 
In fact, the robustness for the effect of high-frequency noise 
exist. With respect to nonlinear double integrator (1), the 
advantage of linear double integrator (5) is its simple imple¬ 
mentation. 

In the next section, the frequency-domain analysis will be 
presented for the nonlinear and linear double integrators. 


III. Erequency-domain analysis based on 

FREQUENCY SWEEP 

In a practical problem, high-frequency noises exist in 
measurement signal a(f). The following analysis concerns 
the robustness behaviors of the nonlinear and linear double 
integrators under high-frequency noises. 

Eor the nonlinear double integrator, an extended version of 
the frequency response method, frequency-sweep method [18, 
19], can be used to approximately analyze and predict the 
nonlinear behaviors of the nonlinear integrator. Even though 
it is only an approximation method, the desirable properties it 
inherits from the frequency response method, and the shortage 
of other, systematic tools for nonlinear integrator analysis, 
make it an indispensable component of the bag of tools of 
practicing control engineers. By frequency-sweep method, we 
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can find that the nonlinear double integrator leads to perform 
precise estimation of integrals and strong rejection of high- 
frequency noise. 

The test of frequency characteristic can be implemented by 
Bode plot htting. For linear or nonlinear double integrators (1) 
or (5), let the input signal be 

a{t) = Ajn sin(a;f) (14) 

where and ui are the amplitude and angular rate of the 
input signal, respectively. Suppose the output of the double 
integrator can be expressed as 


yit) 


Af sin(a;f + Lp) 

\ sin(a;f) cos(wf) 


A f cos ip 
\_ A f sin ip J 


(15) 


where Af,oj and (p are the amplitude, angular rate and phase 
of the output signal, respectively. Let t = 0, h, 2h, • • •, nh, 
and 


2/(0) y{h) ••• y{nh) ] 
sinjwO) sinjwh) ■ • • sin(wn/i) 

cos(a;0) cos{ujh) ■ ■ ■ cos{ujnh) 

-f cos ip, C 2 = Af simp (16) 

where h is the step size. From (15) and (16), it follows that 

y = T-[ci C 2 (17) 

From the least square method, ci and C 2 can be obtained as 
follow: 


= 

4 /^ = 

Cl = 


[ Cl C2 ]^ = (18) 

For the angular rate w, the amplitude and phase of the output 
signal are, respectively, written as 

+ c|, ip = arctan(c 2 /ci) (19) 

Therefore, the amplitude frequency characteristic can be 
described as 

M = 201g(i/M„) = 20\g{^ cl + cl/ A^) (20) 

and the phase frequency characteristic is the phase error 
between the output and input, and it can be described as 


In the following, using frequency-sweep method, we will 
analyze the effects of the observer parameters on the estimate 
performances and robustness. 

A. Frequency characteristics with different e and as 

For the nonlinear and linear double integrators (1) and (5), 
the parameters are selected as follows: ki = 0.1, = 0.1, 

ks = 1; Am = 1; as = a = 0.3,0.5,1, respectively; 
R = 1/e = 3,4,5, respectively. The Bode plots of the 
frequency-domain characteristics with different e and as = a 
are described in Figs. 1(a), 1(b) and 1(c), respectively: Fig. 1(a) 
presents the frequency characteristic of signal tracking; Figs. 
1(b) and 1(c) present the frequency characteristics of onefold 
and double integral estimations, respectively. We can hnd that 
high-frequency noise can be reduced sufficiently. Moreover, 
decreasing parameter e, the cut-off frequency become larger. 
The smaller e is, the signal in wider frequency bandwidth 
can be estimated. However, much noise will pass through the 
observer. On the other hand, increasing parameter e, the cut¬ 
off frequency become smaller, much noise will be rejected. 

Importantly, from Figs. 1(a), 1(b) and 1(c), changing param¬ 
eter as G (0,1], the amplitude frequency characteristics almost 
don’t be affected, but the phase frequency characteristics 
are affected obviously: when as approaches to 1, the phase 
frequency characteristic curves decay slowly near the cut¬ 
off frequency, and phase delay and chattering exist. On the 
other hand, decreasing parameter as G (0,1), the phase 
frequency characteristic curves decay rapidly at the cut-off 
frequency. Relatively smaller as G (0,1) can obtain more 
precise estimations and stronger robustness. 

B. Frequency characteristics with the change of Am 

For nonlinear double integrator (1), the parameters are 
selected as follows: as = 0.3; R = 1/e = 3; fci =0.1, 
k 2 = 0.1, ks = 1; Am = 5,1, 0.5, respectively. The Bode plots 
of the frequency-domain characteristics with the change of Am 
are described in Figs. 2(a), 2(b) and 2(c), respectively: Fig. 2(a) 
presents the frequency characteristic of signal tracking; Figs. 
2(b) and 2(c) present the frequency characteristics of onefold 
and double integral estimations, respectively. It is found that, 
when the magnitude of input signal Am is larger, the cut¬ 
off frequency is relatively smaller, and much noise is reduced 
sufficiently; when the magnitude Am is smaller, the cut-off 
frequency is relatively larger, and the signal in wider frequency 
bandwidth can be estimated. 


'Fe = 'Font - Fin = arctan(c 2 /ci) (21) 

The angular frequency sequence {ook}, where k = 
1, • • • , n, is selected in the interested frequency bandwidth. For 
each angular frequency the frequency bandwidth, the above 
frequency-sweep method is adopted to obtain the values of 
amplitude and phase, respectively. Accordingly, the Bode plots 
of the frequency-domain characteristic can be described. 

In this frequency-domain analysis, oj = 27r/, where / = 
0.1 : 0.5 : 100; h = 0.001; fc = 1 : 1 : 50000. Then input 
signal is a{k) = Am sm{2Trfkh). 


Remark 1: Comparing with ideal integral operators 1/s 
and 1/s^, not only the double integrators can obtain their 
estimations precisely, but also the high-frequency noise is re¬ 
duced sufficiently. Parameter £ affects the low-pass frequency 
bandwidth: Decreasing the perturbation parameter e, the low- 
pass frequency bandwidth is larger, the estimation precision 
becomes better, and relatively higher frequency noise can be 
reduced; on the other hand, increasing perturbation parameter 
£, the low-pass frequency bandwidth is smaller, much noise 
can be reduced sufficiently (See the cases of i? = l/£ = 3,4, 5 
in Figure 1, respectively). Parameter as G (0,1] affects the 
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decay speed of frequency characteristic curves near the cut¬ 
off frequency (See the cases of as = a = 0.3,0.5,1 in 
Figure 1, respectively): Smaller as € (0,1] can obtain more 
precise estimations; Larger as G (0,1] can restrain much 
noise, however, estimation delay phenomena exist. 

IV. Computational analysis and simulations 

In this section, simulation results are presented in order to 
observe the performances of the nonlinear and linear double 
integrators. We consider the simulations of the following case: 
Double integrators for a input signal with high-frequency 
noise. In the simulations, the function of —0.1 x 3.14^ x 
sin(3.14f) is selected as reference signal aos{t). Therefore, 
oo 2 = Jq aos(o')da = 0.1 x 3.14 x cos(3.14f), and oqi = 
fo fo (<t) dadr = 0.1 sin(3.14f). 

Here, the following high-frequency noise 6{t) is selected 
(See the noise in Fig. 3(a)): O.lsin(lOf) + O.lcos(lOf) + 
0.05 sm(50f) -f 0.05 cos(50f). 

Therefore, the input signal a{t) = aos{t) -f 6{t). From the 
frequency-domain analysis, the observer parameters can be 
selected as follows: R = 1/e = 5, fci = 0.1, = 0.1, 

ks = f, as = a = 0.3, a 2 = 2 ^^^^ ~ 3 - 2 a '' initial 
value of the observer is (xi ( 0 ) = 0 ,X 2 ( 0 ) = 1 , 0:3 ( 0 ) = 0 ). 
In the double integrators (1) or (5), xs tracks signal aosit), 
X 2 and xi estimate the onefold and double integrals of signal 
aos{t), respectively. 

Signal 003 (f) tracking, the onefold and double integral 
estimations in 20 seconds are presented in Fig. 3. Fig. 3(a) 
provides signal 093 (f) with stochastic noise. Fig. 3(b) de¬ 
scribes 003 (f) tracking. Figs. 3(c) and 3(d) present the onefold 
and double integral estimations, respectively. Figs. 4(a)-4(d) 
describe the onefold and double integral estimations in 2000 
seconds. 

From the above simulations, despite the existence of the 
intensive high-frequency noise, the nonlinear double integrator 
showed the promising estimation ability and robustness. Fur¬ 
thermore, from Figs. 4(a)-4(d), no drift phenomenon happened 
in the long-time estimations. 

Figs. 5 and 6 describe the estimations of onefold and double 
integrals by the linear double integrator (5) (when 0:3 = 1), in 
20s and 2000s, respectively. From Figs. 5 and 6, the obvious 
estimation delay and slow convergence of double integral exist. 

V. Conclusions 

Based on frequency sweep method, frequency-domain anal¬ 
ysis is proposed for the nonlinear double integrator and 
the linear integrator. All the two types of integrators have 
the strong robustness for the effect of high-frequency noise. 
Comparing to the linear double integrator, the nonlinear in¬ 
tegrator has the better estimation performance and stronger 
robustness. Importantly, the integrator parameters are easily 
to be chosen from the frequency-domain analysis with respect 
to the analysis in time domain. 
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1(a) 



1(b) 



1(c) 


Fig. 1. Frequency-domain characteristics with different e and a. 1(a) Signal 
tracking. 1(b) Onefold integral estimate. 1(c) Double integral estimate. 



2(a) 



2(b) 



2(c) 


Fig. 2. Frequency-domain characteristics with the change of Am - 2(a) Signal 
tracking. 2(b) Onefold integral estimate. 2(c) Double integral estimate. 
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3(a) 
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time{s) 


3(b) 



3(c) 





4(c) 



0 500 1000 1500 2000 

time(s) 


3(d) 


4(d) 


Fig. 3. Integral estimations in 20s (as = 0.3). 3(a) Input signal with noise. 
3(b) Signal tracking. 3(c) Onefold integral estimate. 3(d) Double integral 
estimate. 


Fig. 4. Integral estimations in 2000s (as = 0.3). 4(a) Input signal with noise. 
4(b) Signal tracking. 4(c) Onefold integral estimate. 4(d) Double integral 
estimate. 
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time(s) 


5(a) 




time{s) 



time(s) 


6(a) 



5(b) 


6(b) 





5(c) 


6(c) 






5(d) 


6(d) 


Fig. 5. Integral estimations in 20s (as = 1). 5(a) Input signal with noise. 5(b) 
Signal tracking. 5(c) Onefold integral estimate. 5(d) Double integral estimate. 


Fig. 6. Integral estimations in 2000s (as = 1). 6(a) Input signal with noise. 
6(b) Signal tracking. 6(c) Onefold integral estimate. 6(d) Double integral 
estimate. 





























































































































































